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^ I Abstract. - For binary mixtures of fluids without chemical reactions, but with 

CN I components having different temperatures, the Hamilton principle of least action is 

able to produce the equation of motion for each component and a balance equation of 
the total heat exchange between components. In this nonconservative Gibbs 
dynamical identity connecting the equations of momenta, masses, energy and heat 
^ I exchange allows to deduce the balance equation of energy of the mixture. Due to the 

unknown exchange of heat between components, the number of obtained equations is 
less than the number of field variables. The second law of thermodynamics constrains 
the possible expression of a supplementary constitutive equation closing the system 
of equations. The exchange of energy between components produces an increasing 
rate of entropy and creates a dynamical pressure term associated with the difference 
of temperature between components. This new dynamical pressure term fits with 
the results obtained by classical thermodynamical arguments in [1] and confirms 
^ ■ that the Hamilton principle can afford to obtain the equations of motions for multi- 

^ . temperature mixtures of fluids. 
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1. — Introduction. 

The theory of mixtures considers generally two different kinds of continua: ho- 
mogeneous mixtures (each component occupies the whole mixture volume) and het- 
erogeneous ones (each component occupies only a part of the mixture volume). At 
least four approaches to the construction of two-fluids models are known. 
The flrst one for studying the heterogeneous two-flows is an averaging method (Ishii 
[2]; Nigmatulin [5]). The averaged equations of motion are obtained by applying 
an appropriate averaging operator to the balance laws of mass, energy, etc., valid 
inside each phase [H [5]. A second approach known as Landau method [U [7] was 
used for the construction of a quantum liquid model and was purposed for the ho- 
mogeneous mixtures of fluids [8j. For the total mixture, the method requires the 
balances of mass, momentum, energy, complemented with the Galilean invariance 
principle and the second law of thermodynamics [S]. A third approach is presented 
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in extended thermodynamics; the mixtures are considered as a collection of differ- 
ent media co-existing in the physical space. This approach is done in the context of 
rational thermodynamics [lO] founded on the postulate that each constituent obeys 
the same balance laws as a single fluid fUX [T2l [13]. The thermodynamical pro- 
cesses must verify the second law of thermodynamics and it is possible to purpose 
phenomenological constitutive equations which allow to obtain the structure of con- 
stitutive or production terms (as Pick's law, Fourier's Law etc.) and to close the 
system of equations. 

There exists a different approach based on the Hamilton principle which is used 
for the construction of conservative (non-dissipative) mathematical models of con- 
tinua. The principle was initiated by Lin [14j, Serrin [I5l and many others to obtain 
the governing equations of one component continua [151 [IHl [IZ] and involves an 
Hamilton action. The variations of the Hamilton action are constructed in terms of 
virtual motions of continua which may be defined both in Lagrangian and Eulerian 
coordinates [TH [16] . 

Here, we use variations in the case of fluid mixtures. The variational approach to 
the construction of two-fluid models has been used by many authors (Bedford & 
Drumheller |18]; Berdichevsky [19]; Geurst [20]; Gouin [21]; Gavrilyuk & Gouin 
[22]; Gouin & Ruggeri [23]). 

To study thermodynamical processes by the Hamilton principle, the entropy of the 
total mixture or the entropies of components are added to the field parameters in- 
stead of temperatures. The Lagrangian is the difference between the kinetic energy 
and an internal potential per unit volume depending on the densities, the entropies 
and the relative velocities of the mixture components (and a potential due to exter- 
nal forces). The internal potential per unit volume can be interpreted as a Legendre 
transformation of the internal energy. In this case, it is not necessary to distinguish 
molecular mixtures from heterogeneous fluids when each component occupies only a 
part of the mixture volume [2ll[25]. Consequently, the terms including interaction 
between different components of the mixture do not require constitutive postulates 
difficult to interpret experimentally. They come from the direct knowledge of the 
internal potential per unit volume. 

The assumption of a common temperature for all the components is open to doubt 
for the suspensions of particles [26] as well as in the mixtures of gases in the early 
universe ^7\. By using the Hamilton principle, the existence of several tempera- 
tures (one temperature for each component) must be associated with the existence 
of several entropies (one specific entropy for each component). That will be the aim 
of this paper: the internal potential per unit volume is a function of the densities, 
the entropies and the difference of velocity between components. 
The plan of the article is as follows: 

In section 2, we formulate an extended form of the Hamilton principle of stationary 
action allowing to produce the governing equations of motion for each component 
of a binary mixture. From the invariance of time, we deduce an equation of the 
exchange of heat between components. With two temperatures (one temperature 
for each component), the system of governing equations together with the balance 
of masses is not closed: Hamilton's principle is not able to get a complete set of 
equations when exchange of energies occurs between components. In this case, we 
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need an additive constitutive equation to close the system. 

In section 3, we obtain a Gibbs dynamical identity. This identity allows to obtain 
the equation of energy for the total mixture. The equation of motion for the total 
mixture and the equation of energy are in divergence form; as in [28j, it is not the 
same for the component equations of motion. 

In section 4, we consider the case of weakly dissipative mixtures and introduce an 
average temperature of the total mixture. The average temperature corresponds to 
a local equilibrium different from the real state, with nonequal component temper- 
atures but with the same total internal energy. The total pressure of the mixture in 
the real state is different from the pressure associated with the mixture at the local 
equilibrium [29l [30] . The entropy variation rate of the mixture must be in accor- 
dance with the second law of thermodynamics which implies an additive constitutive 
equation for the pressure. The constitutive equation depends on the physical prop- 
erties of the two constituents and we focus on the fact that the new pressure term 
does not obey the same rule than the other terms due to dissipation. 
In section 5, as an example, we reconsider the special case of a mixture of perfect 
gases. 

Finally, we compare our results with the conclusions obtained in pj by using classical 
arguments of rational thermodynamics. 



2. — Governing equations in conservative cases. 

In a Galilean system of coordinates, the motion of a two-fluid continuum can be 
represented by two diffeomorphisms 

Z„ = $„(z), (a = 1,2) 

or 

X = t and Xq = (t, x), 

where z = (t,x) denotes Eulerian coordinates in a four-dimensional domain u in 
the time-space and = (A, X^) denotes Lagrangian coordinates of the component 
a in a four-dimensional reference space cj^. The conservation of matter for each 
component requires that 

p„ det Fa = p„o (Xa) with Fa = (1) 

where p„o is the reference density in Ua and det ((9x/9Xq) is the Jacobian deter- 
minant of the motion of the component a of density In differentiable cases 
equations ([1]) are equivalent to the equations of density balances 

^ + div(p„v„) = 0, (2) 

where Vq, denotes the velocity of each component a. 
The Lagrangian of the binary system is 

\ /I 

Pa V„ - Pa^a - ^(Pl, P2, ^1, S2, u) 



a=l 
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where the summation is taken over the fluid components (a = 1, 2) and Sa are the 
specific entropies, u = V2 — vi is the relative velocity of components, Qa are 
the external force potentials, r/ is a potential per unit volume of the mixture. The 
Lagrangian L is a function of Pa,^a, Sa and we introduce the quantities 

^«-|^=iv^-#-^., (3) 

(5) 
d L Otj 

where ^ denotes the transposition and 7: — , —- are linear forms. Equation (I5|) 

defines the temperatures (a = 1,2) which are dynamical quantities depending 
on Pi, P2, si, S2 and u. 

To obtain the equations of component motions by means of the Hamilton principle, 
we consider variations of particle motions in the form of surjective mappings, 

Xq, = £iQ,(t, x;Xq,), 

where scalars Xq, are defined in a neighborhood of zero; they are associated with a 
two-parameter family of virtual motions. The real motions correspond to = 
such that Ha(t, x;0) = x); the associated virtual displacements generalize what 
is obtained for a single fluid [HI 



9x, 

The Hamilton action is 



_ dSa{t,^]Ha) 
Oa-^a — 7\ \xa=0 ■ l^j 



L dvdt. 

We first consider the Hamilton principle in the form 

Satt = { - — ) = 6a I L dv dt = 

dXn 



under constraints ([T]), with 6 ad being the variations of a associated with equation 
©. 

From the definition of virtual motions, we obtain in Appendix A the values of 
5q,Vq (x, t), SaPa (x, t) and SaSa (x, t) whcrc (5at;(t,x) is the variation of v at (t, x) 
fixed. By taking into account the formulae in Appendix A and the definitions ([3}[ 

we get 

6aa = j (^RaSaPa + Pa ^a ^a^a - PaTa^aSo^ dv dt 
= I yRa diVc,(p^o ^a'^a) - PaO ^(5aX„) - p^Q Ta K^a\ dVa dt 



A 



where diVa is the divergence operator with respect to the coordinates Xq. In the 
last expression all quantities are considered as functions of [t, Xq,); the functions are 
assumed to be smooth enough in the domain Ua and 6a^a = on its boundary. 
Hence, we get 

and we obtain the equations of component motions in Lagrangian coordinates, 



dy ° ' sx„ "ax. 



a 



where is defined by equation (jl]). 

d F O'v 

By taking into account the identity ° " — = and for A = t, we rewrite 

dt ox 

the equations in Eulerian coordinates. 

The covector is an essential quantity; indeed, p^ko, (and not PqVq, ) is the 
momentum for the component a of the mixture. 

To obtain the equation of energy, we need a second variation of motions associ- 
ated with the time parameter. The variation corresponds to a virtual motion in the 
form 

A = (fit; k), 

where scalar k is defined in a neighborhood of zero. The real motion of the mixture 
corresponds to k = such that (p{t] 0) = t; the associated virtual displacement is 



6X 



dip(t; k) 
Ok 



\k=0 



For a single fluid, the entropy is defined on a reference space Uo associated with 
Lagrangian variables; in conservative motion the specific entropy is conserved along 
the trajectories and in the reference space the entropy depends only on Lagrangian 
variables X and not on A. 

In multi-component fluids, due to exchanges of energy between the components, the 
entropies cannot be conserved along component paths; in the reference spaces cUa, 
the specific entropies Sa depend also on A 

The variation of Hamilton's action associated with the second family of virtual mo- 
tions yields 

6a = 6 I L dxdt = I — 6X dvdt = 0. 
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Q=l 0=1 

d s ds ds 

" ° = — ^ + ^t-^Vq, is the material derivative with respect to velocity Vq,. We 
dt ot ox. 
obtain for the total mixture 

d.s. 



E ^ (y. 



dt 
0=1 

Due to equations ([2]) we obtain the equivalent form 

2 



= with = ^ Q^" + div(p^SQVc,) j T^,. 



(9) 



Equation ([9]) expresses that the exchange of energy between components has a null 
total amount. 



3. — Gibbs dynamical identity and equation of energy. 

Let us prove that equation ([9]) leads to the equation of energy of the mixture. 
We introduce the quantities M^, Ba, S and E such that 

M„ = Pa—jr- + Po^K^ - Pa^r: PaTc.- 



dt ' " 5x dx " ax 



Ba = -g^ + div(p„v„), 
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"S* — ^ ^ Q«, 

0=1 

a=l ^ ^ ^ ^ 

where ri — — u = / is the Legendre transformation of ri with respect to u and 
ou 

corresponds to the volume internal energy of the mixture. We prove in Appendix B 
the following property: 

Theorem: For any motion of the mixture, we have the algebraic identity 

E - Mj v„ + (k^ v„ - i?, + TaSa) B^j-S = 0. (10) 

Relation (fTOll is the general expression of the Gibbs identity in dynamics. Analogous 
identities were obtained earlier for thermocapillary mixtures |2l] and bubbly liquids 
[3T] . Due to the equations of balance of masses ([2]), momenta (JTj) and energy ([9]), 
deduced from Hamilton's principle, which are respectively 

S„ = 0, M„ = (a = 1, 2) and S = 0, 



fi 



we obtain from identity (fTOj) : 

Corollary: The motions of a mixture satisfy the equation of energy balance in the 
form 

^ Wt ^ ^ ^ (^"'^" - Ra) ^ - = • (11) 



Equation ( jlli) appears as the equation of energy when / is the total internal energy. 

The equations of component motions are not written in divergence form. Nev- 
ertheless by summing equations (171) in the form = and taking into account 
equations ([1]) in the form Ba = 0, we obtain by a calculation similar as in [231 [211 [25] 
the balance equation for the total momentum in a divergence form: 



+divfp,v„kl+fp,|!^-r^)l') +p„^ = 0, (12) 
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where I is the identity tensor. In the following, pv = ^^p^k^, = ^^p^v^ is the 



a=l a=l 



total momentum and p = p„ is the mixture density. 



a=l 

System (f[2l).f[71). f[TTl) ). consequence of the Hamilton principle, is a non closed 
system of equations. In a single conservative fluid, the system of motion equations 
is closed by the entropy conservation. In case of mixtures with two entropies, the 
Hamilton principle is not able to close the system of motion equations; we need 
additional arguments to obtain the evolution equations for each entropy Sa by con- 
sidering the behaviors of Qa- 

A possibility to close the system of equations is to consider the case when the mo- 
menta and heat exchanges between the components are rapid enough to have a 
common temperature. This case is connected with a conservative equation for the 
total specific entropy [2lj. 

Another possibility, used by Landau for quantum fluids [7], is to assume that the 
total specific entropy s is convected along the first component trajectory 

dp,s 

+ div (PiS vi) = 0. 

In this case, the constitutive functions are p^^, s, vi, pg, V2, where p^ (a = 1,2) are 
submitted to the constraints and the case of Helium superfluid is a special case of 
our study corresponding to si = s and S2 = 0. Such an hypothesis is not acceptable 
for classical fluids. These assumptions are not valid for heterogeneous mixtures 
where each phase has different pressures and temperatures [26l [27] . 
In the following we consider the case when the mixture is weakly out of equilibrium 
such that the difference of velocities u and the difference of temperatures T2 — Ti 
are small enough with respect to the main field variables. 
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4. — Mixtures weakly out of equilibrium. 

For the sake of simplicity, we neglect the external forces. Generally, the volume 
potential rj is developed in the form [221 ISlj^ 

viPi, P2> Si, S2, u) = e(pi, Si, S2) - b{p^, Si, S2) u^, 

where 6 is a positive function of p^, p2, Si, S2- Properties of convexity of the function rj 
are studied in |25j. When |u| is small enough, the equations of motions are hyperbolic 
|25j . We consider the linear approximation when |u| is small with respect to |vi| and 
|v2|. In linear approximation the volume potential is equal to the volume internal 
energy e, 

V{P1,P2,SI,S2,U) ^ e(pi,P2,Si,S2) = pe{p^,p2,Si,S2), 

where e denotes the internal energy per unit mass. Let us note that the diffusion 
vector j = Pi(vi — v) = P2('^~'^2) is a small momentum vector deduced respectively 
from velocities and densities of the components. The equations of density balances 
can be written in the form 

^ + pdivv = and p^ + divj = 0, (13) 

1 Pi 1 , ■ r ^ , d d d . 

where c = — denotes the concentration 01 component 1 and ^ = tt + tt- ■ v is 
p at ot ax 

the material derivative with respect to the average velocity of the mixture. 

The divergence of a linear operator A is the covector divA such that, for any constant 

vector a, (div A) a = div (A a) and we write Vq,v^ = Vq ® v^,. 

Let us denote by ha = — — the specific enthalpy of the component a. 

For processes with weak diffusion, the equations of component motions get the form, 

d V 

p„ = °^ " + div(p^v„ ® \af = PaTa grad - p^ grad h^. 
The equation of total momentum (fT^ is reduced to 

+ div ( ^ (p^v„ ® Va) - t 

\a=l 

where t = Yla=i is the total stress tensor such that 

2 

tau^ = -Pa Su^, with Pa = PPa^,p^ = Pa(^,P^ ~ ^ ' ^ = X] ^° 



a=l 



The equation of energy (flTl) writes in the simpler form 



d_ 

di 



{^ + ^\ Pa^l^ + div V + ^ Q pavl - V, 



0. 



■^In 1 , the internal energy is the sum of the internal energies of the components 



The internal energy is a natural function of densities and entropies. Due to equation 
©, 

PiTi = p-^{Pi,P2,si,S2) and P2T2 = p —{p^, P2, si, S2). (14) 

Let us denote by e the expression of the specific internal energy as a function of 
p, c, si, S2 such that e{p, c, si, S2) = eipi, P2, si, S2); we get 

de de dp ^ de dc ^ de dsi ^ de ds2 
^ dt ^ dp dt ^ dc dt ^ dsi dt ^ ds2 dt 

d's ds 

Due to the fact that p^ tt — V ^^'^ tt = ^1 ~ ^2? we obtain 

op dc 

de p dp N rr ds2 , . 

pdrt = -pdi^p^^'-^'^drt + p'^'ii+p^^'iF- 

By taking into account that 

and by using equations ([H]), dlSD, equation (ITBi) yields 

de 

p — + p div V + (/ii - h2) div j + (Ti grad Si - T2 grad S2) ^ j = 0. (16) 



Due to equations flMj) . the internal energy can be expressed as a function of 
densities and temperatures of components 

e(Pi> Pi-, Ti, T2) = £:(pi, P2, si, S2). 

As we did in yy, we define the average temperature T associated with Ti and T2 
through the implicit solution of the equation 

e(Pi, P2, T, T) = e(pi, P2, Ti, r2). (17) 

We denote by = — T the difference between component and average tem- 
peratures which are non-equilibrium thermodynamical variables. Near equilibrium, 
equation f|T7|) can be expanded to the first order; then 

^c^e« = with c: = ^{p,,p2,T,T). (18) 



Due to the fact that 

2 



p = ^ p„ Tc, c/sc, + dp^ 

then 

-Qjr{pi,p2,T,T) and pc„=T2^ 

a=l ^ a=l 

9 



P(i = T'^p^-^{pi,P2,T,T) and pel = T^^p^—^ip^, p2,T,T). (19) 

1 1 2 



The definition of tlie total entropy s of the mixture is 

2 

pS = '^PaSa{pl,p2,Tl,T2). (20) 
a=l 

The first order expansion of equation (!20i) yields 

0=1 ^ ^ 

Due to Relations ([HD, ^ 

2 

PS = ^PaSa{pl,p2,T, T) 

and the specific entropy s does not depend on 9i and 62 but only on Pi, P2 and T. 
We denote by e the internal specific energy as a function of p, c, T: 

£{p,c,T) = e{p^,p2,T,T), 

which satisfies the Gibbs equation 

Tds = de — ^dp + (/in — Pi) dc 

where Po (p, c, T) is the equilibrium pressure at temperature T and P2~Pii difference 
of component chemical potentials, is the chemical potential of the whole mixture. 
By taking into account of equation (fT3l) . we get 

de . ^ . ^ ds 

P^+Po dwv + [pi- P2) ^^^J~^^^ = °- 

Moreover, 

pf-Ep«^+d-ife-oj]. (21) 

a=l 

Equation (12T!) yields the relation between the material derivatives of entropy s 
and entropies Si and S2- By taking into account of these results in equation (fT6|) 
and e{p, c, T) = e{pi, P2, Si, S2), we obtain 



d 

T^Pa + (P - Po) div V + ({hi - h2) - (/ii - /ig) + T (ss - si)) div j 



Pa 

a=l 



+ (Gigradsi-esgradsa) j = 0. (22) 

The differences of temperatures Gi = Ti — T and B2 = T2 — T are small with 
respect to T and j is a small diffusion term with respect to the mixture momentum 
pv; consequently, in an approximation to the first order, the term 

(61 grad Si - 62 grad S2) ^ j 

in 



is negligible. Let us consider 

K = ({hi - hi) - (/ii - /ia) + T (s2 - Si) j divj ; 

we get 

K = ((^1 - TiSi) - (/i2 - T2S2) - (fii - ^2) + ©i-si + ©2^2) divj . 

In an approximation to the first order, the term (OiSi + 02^2) divj is negligible. 
Due to the fact that ^^{p^, p2,Ti,T2) = h^ — TaSa is the chemical potential of the 
component a, when j is a small diffusion velocity with respect to average velocity 
V, the term 

(/^i(Pi, P2, Ti,T2)- /i2(Pi, P2, Ti,T2)- (/ii (pi, P2, T, T) - P2 (Pi, P2, T.T))^ divj 

is vanishing in an approximation to the first order. 

Consequently, in an approximation to the first order, equation fl22|) reduces to 

XlPa^^ = (P-Po) divv. (23) 

a=l 

The exchange of energy between components must obey the second law of thermo- 
dynamics: the total entropy rate is an increasing function of time and we consider 
the second law of thermodynamics in the form 

E + div(p„s.v„)] > (24) 

Due to relations ([2]) the Clausius-Duhem inequality flMl) is equivalent to 

2 

V c 
~ dt 



a=l 



Relation (123|) implies that the second member must be positive. Therefore, as usual 
in thermodynamics of irreversible processes, the entropy inequality requires 

IT = p — po = —A div V . (25) 

This expression defines the Lagrange multiplier A of proportionality such that A > 0. 
The dynamical pressure vr is the difference between the pressure in the process out of 
equilibrium with different temperatures for the components and the pressure of the 
mixture assumed in local thermodynamical equilibrium with the common average 
temperature T. Let us notice that equations (p[25!) allow to obtain Qa values. In 
fact, 

p,T(T2-ri)^ = AT2 (divv)^ and p,r(ri - r2)^ = A (divv)^ 
dt dt 

and the system of field equations is closed. 
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5. — Special case of mixture of perfect gases [13j. 

The internal energy of the mixture is the sum of the internal energies of the 
different gas components. We represent these energies as function of density and 
temperature of components 

Hpi,p2,Ti,T2) =pi£i(pi,Ti) + p2e2{p2,T2)- 

Then, 

where pi and p2 are the pressures associated with ei and £2- An expansion to the 
first order in Ti — T = Oi and T2 — T = B2 yields 

p = p^{p„T) + P2{P2,T) + ^y,T) 61 + §r(P2,T) 62. 

From the definition of the average temperature T we obtain as in [29], 

Pl£l(Pl,r) +P2^2(P2>^) = Pl^l(Pl,^l) +P2^2(P2>^2)- 

Consequently, from equation (l25l) . we get 



vr = ^(Pi,T)ei + ^(p2,T)e2. 

Let Ti = T + /3 0, T2 = T + (1 + /5) 9, where 9 = T2 — Ti, an expansion of equation 
(17) to the first order yields the value of j3 

Ci'\p,,T)l3e + Ci'\p2,T){l + l3)e = 0, with = ^(p„, J). 



- a 



Consequently, when pa = kaP^Ta, we obtain 



TT 



PlCv -|- P2 Cv 



In accordance with results obtained by Ruggeri & Simic [30] and Gouin & Ruggeri 
[23], to verify the Clausius-Duhem inequality ( l24l) . 9 must be in the form 



9 = Lt(7i -72)c^it^v with Lt = M ^-^-^ (p^Cf ^ + paC^^^) and M>0 

P2^v 



where 7^, is the ratio of specific heats of component a. 



6. — Conclusion. 

The method by Hamilton can be easily extended to multi-component mixtures 
with multi-temperatures. We obtain the equations of component motions and the 
equation of the total mixture energy. The entropy is not conserved and the sec- 
ond law of thermodynamics reveals the existence of a new dynamical pressure term. 
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As diffusion is a property of ffuid mixtures with different component velocities, the 
dynamical pressure term is a property of fluid mixtures with different component 
temperatures. The dynamical pressure can be measured with the change of volume. 
In the special case of mixture of gases, the dynamical pressure term comes from the 
fact the gases are molecular ly different. 

The Hamilton principle points out that the dynamical pressure can be obtained by 
neglecting viscosity, friction or external heat fluxes. This is a main property of mix- 
tures with multi-temperatures and this fact may have some applications in plasma 
of gases and in the evolution of the early universe |32j . 

In Appendix C, we highlight that constitutive equations for diffusion, viscosity and 
heat flux for mixtures without chemical reaction are consequence of dissipative terms 
whereas the dynamical pressure term can exist with different component tempera- 
tures even if the bulk viscosity is null. 

The results are in complete accordance with the ones by Ruggeri & Simic [20] and 
Gouin & Ruggeri [1]. This is an important veriflcation of the fact that the Hamilton 
principle can be extended to nonconservative mixture motions when components 
have different temperature. A difference with classical thermodynamics methods 
is that the volume internal energy is not necessary the sum of the volume internal 
energies of the components. In this paper, the volume internal energy is a nonsepa- 
rate function of densities and entropies (or temperatures) and is consequently more 
general than in [1] and [30] . 



Appendix A. 

The deflnition of Lagrangian coordinates X^, implies " + ^ " Vq, = 0. By 

at ox 

taking the derivative with respect to Xa, we obtain the following equation for virtual 
displacements (equation ([6])) associated with the flrst virtual motion family 

^ ^ K + da^a = 0. 

ot ox ax 

Then, we get 



5aVa (X, t) = - (5aXc 



Equation ([T]) yields 



dPao 



6^p^ (x, t) det (x, t) + (det F,) = <5,X,. (26) 
By using the Euler-Jacobi identity 

(5„(det F,) = det F,(x,t) tr (^F'' 5^F^ 

with 

So^F^ (x,t) = - F„ (x,t) (x,t) (x,t) 
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and 

we deduce 



5«(det F„ ) = - det F„ tr (^So^P-' F^^ = - det tr(^^^y 



or, 



5Q,(detFQ,) = — det Fq, diVa{Sa^a)- (27) 
By substituting equation (127|) into equation (12^ . we obtain 



CA.Q, 

Appendix B. 

The proof of the Gibbs identity is obtained by summing the following algebraic identities 

a — e, 

For the external potentials 0,a, 
a. 

— h dZV [pjla^a) - Pa " BaUa " Pa = 

For the velocity fields Vq, 
b. 

^ Q Pa + (^Pa V„ (^vl 

Let us introduce i = — — — . Then 

on 

\ / a=l ^ 

and the three following identities c — e prove the formula 
c. 



dvi dp ( dfj \ d f Otj \ Otj [dp 

" +div( -— Pa'^a] - Pa TT " 1^ ^ + (Pa ^a) ) = 0, 



dp^ dt \dp^ 5x \dp^J " V 5* 

d. 

m (^^a I rp dSa rn ^otSa ^ 

P^Ta^ + Pa T-^^- - P^T'^^ = 0' 



u 



e. 



dt 



+ div 



a=l 



ft.|((-ir-)+..(-ir-^iv„ 

dt \ ) p„ 5x 



- (-11 



dt 



+ div Va) = 0. 



Appendix C. 

We consider a more general case of a mixture when the Hamilton principle cannot be 
applied. This case consists of a weak dissipative process with diffusion, viscosity and heat 
transfers. The balance of masses, momenta and energy are simply expressed by adding 
dissipative terms to the expressions obtained in section 4 



Ba = 0, Mi = {a = 1, 2) and E'^ = 0, (28) 



such that 



Ba = ^ + div(p,v,). 



= - (div cr^)^ - 

2 

^'^ = E + ^divq«-cr^v„. 



Q = l 



On the right hand side, qo, is the heat flux vector, nia is the momentum production and 
cr^ is the viscous part of the stress tensor of constituent a. Due to the total conservation 
of momentum of the mixture, Yl'a=i = [T2]. 
The dynamics Gibbs identity (|lUp can be transformed as 

with 

2 

S'^ = 5 + divq+^m>,-tr (o"^ D„), 



a=l 



where Da = ^ I + 



^) andq= ^q«. 

^ / a=l 



9x \ 9> 

The second law of thermodynamics is expressed in the form 

i:^«^ + div|^>0. (29) 

In the second order approximation, with small external heat fluxes qo- and small difference 
of temperature Ti — T2, equation ([29|) is equivalent to 

2 , 

P«^ + div->0. 
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If we write m = mi = 



— m2, we obtain by calculations similar to section 4 



Pr. —r— + div - 



.a=l 



dt 



+ s = o, 



where 



S = (p — Po) div ^r +^ grad T — m^'^u — ^ tr (cr^ D^) < 



is the entropy production. 

Classical methods of thermodynamics of irreversible process ( TIP) yield equation (j25p for 
the dynamical pressure term together with Fourier and Navier-Stokes laws [11] . 
Term m^u yields the coefficient x of proportionality such that x ^ ^-iid 



m 



-Xu = x(vi - V2). 



For slow isothermal motions, the difference between the components of equations ([281 
yields in an approximation to the first order 



Ml M2 

Pi P2 



grad(^i - ^^2)■ 



Here 



drj 
dpc 



TaSa denotes the chemical potential of component a at temperature T. 



By neglecting the viscous terms, we obtain 
P 



grad/U 



P1P2 



m [p = Pi — P2) or grad/i = — KU with k 



P1P2 



■X, 



which is an expression of the Fick law. 

Therefore, in this formulation, our results coincide with the ones obtained by arguments 
of classical thermodynamics [1]. 
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